Normal forms and entanglement measures for multipartite quantum states 
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A general mathematical framework is presented to describe local equivalence classes of multipar- 
tite quantum states under the action of local unitary and local filtering operations. This yields 
multipartite generalizations of the singular value decomposition. The analysis naturally leads to the 
introduction of entanglement measures quantifying the multipartite entanglement (as generalizations 
of the concurrence and the 3-tangle), and the optimal local filtering operations maximizing these 
entanglement monotones are obtained. Moreover a natural extension of the definition of GHZ-states 
to e.g. 2 x 2 x TV systems is obtained. 



One of the major challenges in the field of quantum 
information theory is to get a deep understanding of 
how local operations assisted by classical communication 
(LOCC) performed on a multipartite quantum system 
can affect the entanglement between the spatially sepa- 
rated systems. In this paper we investigate this problem 
in the case that only operations on one copy of the sys- 
tem are allowed. This is different from the general setup 
of entanglement distillation, where global operations on a 
large (infinite) number of copies are performed to concen- 
trate the entanglement in a few copies. The main motiva- 
tion of this work was to characterize the optimal filtering 
operations (SLOCC) to be performed on one copy of a 
multipartite system such that, with a non zero chance, a 
state with maximal possible entanglement is obtained. In 
other words, we want to design the optimal filtering oper- 
ations for a given state, such that with a certain chance 
we prepare the optimal attainable one. Of course this 
leads to the introduction of local equivalence classes. 

In the case of a pure state of two qubits, this opti- 
mal filtering procedure is commonly known as the Pro- 
crustean method Q. Following the work of Gisin Q, 
Horodecki 0] , Linden et al. 3 and Kent et al. 0, E| , the 
optimal filtering procedure for mixed states of two qubits 
was recently derived in • In this paper we extend these 
ideas to pure and mixed multipartite systems of qudits 
of arbitrary dimension. 

The optimal filtering operations in Q were derived by 
proving the existence of a decomposition of a mixed state 
of two qubits as a unique Bell diagonal state multiplied 
left and right by a tensor product representing local op- 
erations. A Bell diagonal state is special in the sense 
that one party alone cannot acquire any information at 
all about the state: its local density operator is equal to 
the identity. This can readily be generalized to multipar- 
tite systems of arbitrary dimensions, and the existence of 
local operations transforming a generic state to a unique 
state with all local density operators equal to the identity 
will be proved. In the case of pure states, this decomposi- 
tion leads to a transparent method of deriving essentially 
different states such as GHZ- and W-states M. 



We then proceed to show that all quantities exhibiting 
some kind of invariance under the considered SLOCC op- 
erations are entanglement monotones |9j ■ It is shown that 
the concurrence and the 3-tangle, introduced by Woot- 
ters et al. belong to this class of entanglement 

measures. Therefore a natural generalization of these 
measures is obtained to systems of arbitrary dimensions 
and an arbitrary number of parties. 

A subsequent part of the paper is concerned with find- 
ing the optimal filtering operations for a given multipar- 
tite state. It is shown that the SLOCC operations bring- 
ing a state into its unique normal form maximize all the 
introduced entanglement monotones. This was expected 
in the light of the work by Nielsen about majorization 
[l2| : the notion of local disorder is intimately connected 
to the existence of entanglement. 

Finally, the appendix presents some results on the 
characterization of local unitary equivalence classes, 
yielding a natural and constructive but non-unique gener- 
alization of the singular value decomposition to the mul- 
tilinear setting. 



Normal forms under SLOCC operations 

Let us first consider the case of pure states. The main 
goal is to study equivalence classes under general local 
transformations of the kind = Ax O • • • ® A n \i/j) with 
{Ai} arbitrary matrices. These kind of transformations 
are called SLOCC transformations 8] (from stochastic lo- 
cal operations assisted by classical communication), and 
are also called local filtering operations. It will turn out 
very useful to restrict ourselves to SLOCC transforma- 
tions where all {Ai} are full rank (remark that entangle- 
ment is lost whenever a Ai is not full rank). For conve- 
nience, we will consider all {Ai} to belong to SL(n,C), 
the group of square complex matrices having determinant 
equal to 1, and consider unnormalized states. 

Let us formulate the following central theorem: 

Theorem 1 Consider an Ni x N2 x • • • N p pure multipar- 
tite state (or tensor). Then this state (tensor) can con- 
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structively be transformed into a normal form by determi- 
nant 1 SLOCC operations. The local density operators of 
the normal form are all proportional to the identity, and 
the normal form is unique up to local unitary transfor- 
mations. Moreover, the state connected to the original 
one by determinant 1 SLOCC operations with the mini- 
mal possible norm (i.e. trace of the unnormalized density 
operator) is in normal form. 

Proof: We will give a constructive proof of this the- 
orem that can directly be translated into matlab code. 
The idea is that the local determinant 1 operators A.- L 
bringing tp into its normal form can be iteratively de- 
termined by a procedure where at each step the trace 
of = p is minimized by a local filtering opera- 

tion of one party. Consider therefore the partial trace 
pi = Tt2 - p (p)- If Pi is full rank, there exists an operator 
X with determinant 1 such that p' x — Xp\X^ ~ In, . In- 
deed, X = IdetGoi)! 1 / 2 ^^)- 1 does the job Ha, and 
we have p[ = def^i) 1 ^ 1 /^. We also have the relation: 

Tr(p') = A r idet(p 1 ) 1 / Ari <Tr( Pl ), (1) 

where p' = (X <g I ■ ■ ■ <g (ip\(X (g /• • • ® J)L This 
inequality follows from the fact that the geometric mean 
is always smaller than the arithmetic mean, with equality 
iff pi is proportional to the identity. Therefore the trace 
of p decreases after this operation. We can now repeat 
this procedure with the other parties, and then repeat 
everything iteratively over and over again. After each 
iteration, the trace of p will decrease unless all partial 
traces are equal to the identity. Because the trace of 
a positive definite operator is bounded from below, we 
know that the decrements become arbitrarily small and 
following equation this implies that all partial traces 
converge to operators arbitrarily close to the identity. 
We still have to consider the case where we encounter a 
Pi that is not full rank. Then there exists a series of X 
whose norm tends to infinity but has determinant 1 such 
that XpiX^ = 0, leading to a normal form identical to 
zero, clearly the positive operator with minimal possible 
trace. This ends the proof of the existence of the normal 
form. 

Consider now the a state that is normal form; then due 
to the construction of the proof, the trace can always be 
decreased by determinant 1 SLOCC operations, unless 
the state is in normal form. 

As pointed out by Briand, Luque and Thibon ^3 > the 
normal form is unique up to local unitaries: the Kempf- 
Ness criterion proves the uniqueness in the case of a 
closed orbit, and there is always a unique closed orbit 
in the closure of an arbitrary orbit [14j. This ends the 
proof. □ 

Let us now return to the general theorem ^ This 
theorem is very fundamental in that it states that each 
pure multipartite state can be transformed into a unique 



state with the property that all local density operators 
are proportional to the identity. States in normal form 
are clearly expected to be maximally entangled states. 
As we will argue later, the normal form is the state with 
the maximal amount of entanglement that can be created 
locally and probabilistically from the original state. 

Let us next prove that the normal form is continu- 
ous with respect to perturbations of the entries of the 
original density matrix p. First of all note that the non- 
uniqueness due to the local unitaries can be circumvented 
by imposing all Ai to be hermitian. The following lemma 
shows that the normal form is robust against perturba- 
tions or noise: 

Lemma 1 // the SLOCC operations bringing the state 
into the normal form introduced in theorem^are chosen 
to be hermitian, and if they turn out to be finite, then the 
normal form is continuous with respect to the entries of 
the state. 

Proof: Let us consider p = (Ai (g • • • g> A p )a(A 1 (g • • • <g 
A p Y and a perturbation p resulting in {Ai} and a. The 
following formula is readily verified: 

(Ai <g • • • <g Ap^piAi <g) • • • (g) A p )^ = 
v 

& + ^2 (( 1 ® • • " A i lA * ■ ■ ■ ® + h - c -) ■ 

i=l 

As all {Ai} are hermitian and have determinant 1, all 
A~ l Ai are skew- hermitian and the second term lives in 
another subspace S% than the first term a who lives in 
subspace S±. & can therefore be obtained by projecting 
the left hand side parallel to onto Si. As p is finite 
and all {Ai} have determinant one and are finite, this 
projection is of course also finite. This proofs that a 
is of the same order of magnitude as p, which ends the 
proof. □ 

Note that we have also proven continuity with respect 
to mixing. 

Let us now discuss some peculiarities. The fact that 
the algorithm can converge to zero despite the fact that 
all Ai have determinant equal to 1 is a consequence of the 
fact that SL(n,C) is not compact: there exist states that 
can only be brought into their respective normal form by 
infinite transformations, although the class of states with 
this property is clearly of measure zero. As an example 
consider the W-state Q = |001) + |010) + |100). The 
following identity is easily checked: 

The normal form corresponding to the W^-state is there- 
fore equal to zero, clearly the state with the minimal pos- 
sible trace. This is interesting, as it will be shown that 
a normal form is zero iff a whole class of entanglement 
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monotones is equal to zero. Therefore the states with 
normal form equal to zero are fundamentally different 
from those with finite normal form, and this leads to the 
generalization of the W-class to arbitrary dimensions. 

It thus happens that some states have normal form 
equal to 0. This also happens if the state does not have 
full support on the Hilbert space in that one partial trace 
Pi is rank deficient. Note that states which do not have 
full support on the Hilbert space, such as pure states from 
which one party is fully separable, all have normal form 
equal to zero. It will indeed turn out that the amount of 
multipartite entanglement present in a state can be quan- 
tified by the trace of the obtained normal form, which is 
clearly zero in the case of separable states. On the other 
hand, the only normalized states that are already in nor- 
mal form are precisely the maximally entangled states: 
in the case of three qubits for example, the only state 
with the property that all its local density operators are 
proportional to the identity is the GHZ-state. 

As a last remark, we give an example of a state that is 
brought into a non-zero normal form by SLOCC opera- 
tors that are unbounded: 

\ip) ~ a(|0000) + |1111>) + |01)(|10) + |01» (2) 

The normal form is just given by the GHZ-state ( 1 0000) + 
1 1111)), but as can be derived from the results presented 
in [l5j, infinite SLOCC transformations are needed to 
reach this. 



Entanglement monotones 

Until now we contented ourselves to characterize the 
orbits generated by local unitary or SLOCC operations, 
but we have not tried to quantify the entanglement 
present in a state. The SLOCC normal form introduced 
in the previous section however gives us a strong hint of 
how to do this. Note that all separable states have a nor- 
mal form equal to zero, and that the known maximally 
entangled states such as Bell-states and GHZ-states are 
the only ones of their class that are in normal form. 

This suggests a very general way of constructing en- 
tanglement monotones: 

Theorem 2 Consider a linearly homogeneous positive 
function of a pure (unnormalized) state M{p = \tp)(ip\) 
that remains invariant under determinant 1 SLOCC op- 
erations. Then M{\ij))(ip\) is an entanglement monotone. 

Proof: A quantity M{p) is an entanglement monotone iff 
its expected value does not increase under the action of 
every local operation. It is therefore sufficient to show 
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{Ax I)p{Ax 2)t 



that for every local A\ < In 1 , A\ = y In x 
holds that 

M(p) > Tr {{Ax 2)p(Ai I) t ) 



A\Ai, it 



,Tr((Ai®I)p{Ai®I)1), 
+Tr {{Ax I)p{Ax ® iY) . 

( {a x ® i)p{Ax ® iy 



M 



Tr {{Ax ® I)p{Ax ® J)t 



If Ax is full rank, it can be transformed to a determi- 
nant 1 matrix by dividing it by det{Ax) 1 ^ Nl . Due to the 
homogeneity of M{ap) = aM{p) the previous inequality 
is equivalent to 



M{p) > (Idet^i)!^ 1 + \det{Ax)\ 2/1 ^)M{p). 

As the arithmetic mean always exceeds the geometric 
mean, this inequality is always satisfied. This argument 
can be easily completed to the cases where A, is not full 
rank due to continuity. The same argument can then be 
repeated for the other Ai, which ends the proof. □ 

Entanglement monotones of the above class can read- 
ily be constructed using the completely antisymmetric 
tensor e^...^. 

Indeed, it holds that A^^A^ ■ ■ ■ A iN j If ej 1 ...j N = 
det(A)ej 1 ...j Jv ., and as det(A) = 1 this leads to invari- 
ant quantities under determinant 1 SLOCC operations. 
These quantities seem to be related to hyperdeterminants 
[l^ ? ] , and those latter seem to be a subclass of the quan- 
tities considered here. 

Consider for example the case of two qubits. The quan- 
tity 
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is clearly of the considered class, and it happens to be 
the celebrated concurrence entanglement measure |llj . 
In the case of three qubits, the simplest non-trivial homo- 
geneous quantity invariant under determinant 1 SLOCC 
operations is given by 

IV'ujlfel "^iijlk^i^h^s, r 4 ) i4j4k 4 (-i 1 i2 e i3ii e ji 02 e j3j4 e k 1 k 3 ^k 2 k i | ^ 

(Note that we use the Einstein summation convention.) 
This happens to the square root of the 3-tangle intro- 
duced by Wootters et al.^fji which quantifies the true 
tripartite entanglement. 

More generally, as the considered entanglement mono- 
tones are invariant under the determinant 1 SLOCC op- 
erations, the number of independent entanglement mono- 
tones is equal to the degrees of freedom of the normal 
form obtained in the case of a pure state minus the de- 
grees of freedom induced by the local unitary operations. 
Indeed, this is the amount of invariants of the whole class 
of states connected by SLOCC operations. It is then eas- 
ily proven that a normal form is equal to zero if and only 
if all the considered entanglement monotones are equal 
to zero: the entanglement monotones are homogeneous 
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functions of the normal form, and if the normal form is 
not equal to zero there always exists an SLOCC invariant 
quantity that is different from zero. 

In the case of 4 qubits for example, parameter count- 
ing leads to (2 • 2 4 — 2) — 4 • 6 = 6 (a state has 32 degrees 
of freedom —2 to an irrelevant phase and the 4 SX(2,C) 
matrices have each 6 degrees of freedom) independent en- 
tanglement monotones. The simplest monotone is given 
by 

14*11 jik 1 l 1 lf}i 2 j2k2h e ili2 e jlj2 e k 1 k2 e lil2 I J {") 

and the other 5 entanglement monotones can be obtained 
by including more factors; an example is 

V^IV'il jlfcl/l 4 , i2j2k2h'4'i3j3k3h4 } i4.j4.k4.U ■ 

e ili2 e i3i4. e hl2 e l3U e jlj3 e j2j4 e kik 3 £k2k4,\ ^ (4) 

These are clearly generalizations of the concurrence and 
the 3-tangle to four parties. Note however that the situa- 
tion here is more complicated due to the existence of mul- 
tiple independent entanglement monotones. Note also 
that there exist biseparable states that can be brought 
into a non-zero normal form by determinant 1 SLOCC 
operations. Consider for example the tensor product 
of two Bell states; all local density operators are pro- 
portional to the identity, the value of the entangle- 
ment monotones (01 and is respectively given by 1 
and 1/V2 (as opposed to 1 and 1 for the GHZ-state 
(| 0000) + |1111))/V2), and nevertheless no true 4-partite 
entanglement is present. 

If the subsystems happen to be of unequal dimension, 
then the respective subdimensions should be chosen not 
larger than the maximal allowed dimension such that all 
local density operators remain full rank. In a 2 x 2 x N 
system for example, a pure state can only have full sup- 
port on the 2x2x4 subspace, and therefore it makes 
no sense to calculate the normal form with TV > 4: one 
can always first rotate the TV-dimensional system into a 
4-dimensional one by local unitary operations, and pro- 
ceed by calculating the normal form for the 2x2x4 
system. More generally, if the dimension of the largest 
subsystem does not exceed the product of all the other 
ones, then generically the normal form will not be equal 
to zero, leading to non-trivial entanglement monotones. 
As an example, consider a 2 x 2 x 4 system; there are 
more local SLOCC parameters than the number of de- 
grees of freedom, so there will be only one entanglement 
monotone (as is the case in the 2x2 and 2x2x2 case). 
The 2x2x4 tangle is given by: 

"v/4/3 1 ^ ] '4'iijik l '4 , i2j2k 2 '4'i3hk3' l Piijik4,- 

1 1/2 

e ili 2 e i 3 ii e jlj3 e j2ji e k 1 k2k3k 4 , \ 

The factor \/4/3 is included to ensure that the state in 
normal form 

(|000) + |011> + |102) + |H3>)/2 (5) 



has tangle given by 1. Indeed, as will be shown in the 
following section, the maximal value of the tangle is al- 
ways obtained for states in normal form, and this is the 
unique state (up to LU) having all its local density oper- 
ators proportional to the identity. Note that this state is 
therefore the generalization of the GHZ state to 2 x 2 x 4 
systems. 

For completeness, let us also give a formula for the 
2x2x3 tangle: 

I ^ ] V'il j'l,fcl "0*2 J2.fc2 4 , i3,j3,k 3 '4'ii ,ji,ki V'is jBjfcsV'isAifo • 

e ilii e i 2 i5 e i3i6 e jlji e 32j5 e j3j6 e kik2k3('k4k !! ,k e | 

The state maximizing this entanglement monotone (the 
number is bounded by 1) is the generalization of the 
GHZ to the 2 x 2 x 3 case: 

_L|000) + -^|011) + -^|101) + -^|112). (6) 

Let us finally give a non-trivial example of an entan- 
glement monotone of the considered class in the case of 
three qutrits: 

iijiki 4 ) %232k2 ^3 j3&3 'tpiijiki Ipinjn fe 5 V'ieie fe 6 
e ili2i3 e i4i5i6 e jlj 2 ji e j3j536 e kik 5 k 6 ^k2k 3 k i | ^ ■ 

The other (2-3 3 — 1) — (3-16)— 1 = 4 independent entangle- 
ment monotones can again be constructed by including 
more factors. 

Optimal Filtering 

A natural question now arises: characterize the opti- 
mal SLOCC operations to be performed on one copy of a 
multipartite system such that, with a non zero chance, a 
state with maximal possible multipartite entanglement is 
obtained. This question is of importance for experimen- 
talists as in general they are not able to perform joint op- 
erations on multiple copies of the system. Therefore the 
procedure outlined here often represents the best entan- 
glement distillation procedure that is practically achiev- 
able. 

In the previous section a whole class of entanglement 
monotones that measures the amount of multipartite en- 
tanglement were introduced. The following theorem can 
easily be proved using the techniques of theorem ^ 

Theorem 3 Consider a pure multipartite state, then the 
local filtering operations that maximize all entanglement 
monotones introduced in theorem\^are represented by op- 
erators proportional to the determinant 1 SLOCC opera- 
tions that transform the state into its normal form. 

Proof: The proof of this theorem is surprisingly simple. 
Indeed, all the quantities introduced in theorem [21 are 




5 



invariant under determinant 1 SLOCC operations if the 
states do not get normalized. The value of an entan- 
glement monotone however only makes sense if defined 
on normalized states, and due to the linear homogeneity 
of the entanglement monotones, the following identity 
holds: 

M ( {^A^p^Arf \ = M(p) 

\Tt ((OiAi)p(®i^)t) J Tr ((® i A i )p(® i A i y) 

The optimal filtering operators are then obtained by the 
{Ai} minimizing 

Tr ((® 4 A l )p(® l A i ) t ) . (7) 

But this problem was solved in theorem ^ where it was 
proved that the {Ai} bringing the state into its normal 
form minimize this trace. □ 

It is therefore proved that the (reversible) procedure of 
washing out the local correlations maximizes the multi- 
partite entanglement as measured by the generalization 
of the tangle. This is in complete accordance with the 
results of majorization |12j . where it is shown that the 
notion of local disorder is intimately connected to the 
amount of entanglement present. Therefore we have sup- 
porting evidence to call pure states in normal form max- 
imally entangled with relation to their SLOCC orbit. 

The mixed state case. 

The normal form derived in theorem ^ can readily be 
generalized to the case where the state is mixed, i.e. the 
case where the density operator is a convex sum of pure 
states. Indeed, nowhere in the proof of the theorem it 
was used that the state p was pure; the same holds for 
the continuity for the normal form. We have therefore 
proven: 

Theorem 4 Consider an N\ x N2 x • • ■ N m mixed mul- 
tipartite state. Then this state can be brought into a nor- 
mal form by determinant 1 SLOCC operations, where the 
normal form has all local density operators proportional 
to the identity, and the normal form is unique up to lo- 
cal unitary operations. Moreover the trace of the normal 
form is the minimal one that can be obtained by determi- 
nant 1 SLOCC operations. If the SLOCC operations are 
chosen to be hermitian, then the normal form is contin- 
uous with respect to perturbations of the original state. 

Note that if p is full rank, its normal form will never 
converge to zero: the determinant of the density operator 
is constant under SLOCC operations. 

It is also possible to adopt the results about entangle- 
ment monotones. First of all we extend the definition 
of an entanglement monotone p p that is defined on pure 



states and that is linearly homogeneous in p by the con- 
vex roof formalism: 

Mm(p) = mm y^PiUpQipi))- (8) 
1 

Here the optimization has to be done over all pure state 
decompositions of the state. The fact that the pure state 
entanglement monotone is linearly homogeneous in p en- 
sures that p m is, on average, not increasing under local 
operations, and therefore assures that p m is an entan- 
glement monotone. Moreover, it is obvious that these 
entanglement monotones are again invariant under de- 
terminant 1 SLOCC operations. The results on optimal 
filtering for mixed states also readily apply, and therefore 
we arrive at the following very powerful result: 

Theorem 5 The local filtering operations bringing a 
mixed state into its normal form are exactly the ones 
that maximize the entanglement monotones that remain 
invariant under determinant 1 SLOCC operations. 

This result is remarkable, because there does typically 
not exist a way of actually calculating the value of an 
entanglement monotone defined on a mixed state: find- 
ing the optimal pure state decomposition of a state with 
relation to the convex roof formalism for a given EM is 
excessively difficult and has until now only been proven 
possible for the concurrence (i.e. the case of two qubits). 
So although we cannot calculate the entanglement mono- 
tone, we know how to maximize it! This particularly ap- 
plies to mixed states of three qubits: we have proven how 
to maximize the 3-tangle, although we don't know how 
to calculate it. 

Note that this optimal filtering procedure produces 
non-trivial results even in the case of two qubits: it proves 
that the concurrence and therefore the entanglement of 
formation of a mixed state of two qubits is maximized by 
the SLOCC operations bringing the state into its unique 
(Bell-diagonal) normal form. 

Conclusion 

In conclusion, we presented a constructive way of 
bringing a single copy of a quantum state into normal 
form under local filtering operations. This normal form 
is such that all local information is washed out (i.e. the 
local density operators are maximally mixed). We pre- 
sented qualitative and quantitative arguments why the 
amount of entanglement of states in normal form can- 
not be enlarged by local operations, and introduced a 
whole class of entanglement measures that are a direct 
generalization of concurrence and 3-tangle to systems of 
arbitrary dimension. This sheds some new light on the 
difficulty encountered in classifying, understanding and 
unravelling the mysteries of multipartite quantum entan- 
glement. 
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APPENDIX: NORMAL FORMS UNDER LOCAL 
UNITARY OPERATIONS 

Consider a general multipartite state with m parties 
defined on a n\ ® n 2 ■ ■ ■ um dimensional Hubert space: 



(9) 



In this appendix, we try to solve the following natural 
question: is there a method to verify if two states \ipi) 
and 1^2) are equivalent up to local unitary transforma- 
tions? In the bipartite case, this problem can readily be 
solved using the singular value decomposition, and we 
therefore ask for some kind of generalization of this diag- 
onal normal form. Let us state the following theorem (see 
also Carteret et al. 17]), which is a weak generalization 
of the SVD: 

Theorem 6 Given a general complex tensor V'ii— i m with 
dimensions ni = ni = • • • = n m = n, then there exist 
local unitaries Ui such that all the following entries in 
the tensor ip' = U\ ® • • • ® U m ipi 1 ...i m are set equal to 
zero: 



VI < 3 < n,Vk > j : ip'i 



. = 



v', , .... = 



3,3,— ,3,j,k 
3,J,---3,k,3 



^3,k,j, 
^,3,- 



,3,3 



= 



,3,3 



= 0. 



Moreover all entries ip' n n ... n i n ... n , i < n can be made 
real and positive. If the number of parties exceeds 2, then 
the normal form is typically not unique up to permuta- 
tions, but there exist a discrete number of different nor- 
mal forms with the aforementioned property. The num- 
ber of zeros however can generically not be increased by 
further local unitary operations. 

Proof: unlike the proof in |17j . this proof is construc- 
tive and can readily be translated into matlab code to 
calculate the normal form numerically. Consider first all 
entries with at least m — 1 times 1 in its indices, and 
define the vectors x\ = ^1,1,1, •••,!) x 1 = ' i Pi,i,i,---,ii ■■■ , 
x™ = V-'i,i.- -,i.i- Define now a recursive algorithm that 
goes as follows: rotate x 1 to Hx 1 !! [1, 0, • • • 0] by a unitary 
transformation, apply the same transformtion on the full 
tensor, and define x 2 — ^1,1,1,— ,1 with ip the transformed 
tensor. Now do the same thing with x 2 , ... x m and then 
again with x , until the algorithm converges. This al- 
gorithm will certainly converge because at each step the 



(1, 1, ■ • ■ 1) entry of ip becomes larger and larger, unless all 
entries (1, 1, • • • , 1, i, 1, • • • 1) are equal to zero; moreover 
its value is bounded above because the unitary group is 
compact. Next exactly the same algorithm can be ap- 
plied to the subtensor of ip defined as the one with all 
entries larger or equal to 2 (it is easy to check that the 
zeros obtained in the first step will remain zero by this 
kind of action). Next we can again do the same thing 
of another (smaller) subtensor, proving that indeed all 
zeros quoted in the theorem can be made. 
It is straightforward to prove that the entries 

di-nin»«i' — n can a ^ ^ e ma de real and positive 
by further diagonal unitary transformations. 
Let us finally prove that no more zeros can be made by 
whatever unitaries (in the generic case). This follows 
from the fact that a unitary n x n matrix has n 2 con- 



tinuous real degrees of freedom, but that only n 



of 



them can be used to produce zeros as the other n degrees 
of freedom can be imbedded in a diagonal unitary with 
just phases. Counting of the number of zeros produced 
indeed leads to 



m m — 1 



2^ max ( n 
3=1 fc=i 



k,o) 



n(n — 1) 



(10) 



which indeed corresponds to the m(n 2 — n) degrees of 
freedom as the zeros are " complex" . 
The non-uniqueness of the normal form obtained is sur- 
prising but can readily be verified by implementing the 
algorithm on a generic tensor; typically the algorithm 
converges to one out of a finite number of possible differ- 
ent normal forms. □ 

As a first example, consider a system of three qubits. 
Unfolding the 2x2x2 tensor in two 2x2 matrices, the 
following entries can always be made equal to zero: 



x 
x 



x 

X X 



(11) 



Here x is used to denote a non-zero entry. In this case, 
it is easy to see that 4 of the remaining 5 entries can be 
made real by multiplying with appropriate diagonal local 
unitaries. This is equivalent to the normal form obtained 
by Acin et al.flij]. 

A more sophisticated example is the 3x3x3 case, 
whose normal form looks like 



(12) 



It is also straightforward to generalize the previous the- 
orem (and constructive proof) to s yste ms with different 
subdimensions (see Carteret et al.|ljj for an existence 
proof); the algorithm of the previous proof can readily 
be extended to this case. Let us for example consider the 




7 



normal form of the N x 2 x 2 case 



fx 






a; 
a; 
a; 







fO x\ 

X X 

x 







(13) 



V \0 Oj \Q Oj J 

This case is of particular interest as it is describes a state 
of two qubits entangled with the rest of the world. 
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